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Abstract 

Introduction: The response to exogenous epinephrine (Ep) is difficult to predict given the multitude of factors 
involved such as broad pharmacokinetic and pharmacodynamic between-subject variabilities, which may be more 
pronounced in children. We investigated the pharmacokinetics and pharmacodynamics of Ep, co-administered with 
milrinone, in children who underwent open heart surgical repair for congenital defects following cardiopulmonary 
bypass, including associated variability factors. 

Methods: Thirty-nine children with a high risk of low cardiac output syndrome were prospectively enrolled. Ep 
pharmacokinetics, hemodynamic and metabolic effects were analyzed using the non-linear mixed effects modeling 
software MONOLIX. According to the final model, an Ep dosing simulation was suggested. 

Results: Ep dosing infusions ranged from 0.01 to 0.23 ug.kg^.mirt 1 in children whose weight ranged from 2.5 to 
58 kg. A one-compartment open model with linear elimination adequately described the Ep concentration-time 
courses. Bodyweight (BW) was the main covariate influencing clearance (CL) and endogenous Ep production rate 
(qO) via an allometric relationship: CL(BWi) = 9 CL x (BWi) 3/4 and qO(BWi) = 0 qO x (BWi ) 3/4 . The increase in heart rate 
(HR) and mean arterial pressure (MAP) as a function of Ep concentration were well described using an Emax model. 
The effect of age was significant on HR and MAP basal level parameters. Assuming that Ep stimulated the production 
rate of plasma glucose, the increases in plasma glucose and lactate levels were well described by turnover models 
without any significant effect of age, BW or exogenous glucose supply. 

Conclusions: According to this population analysis, the developmental effects of BW and age explained a part of the 
pharmacokinetic and pharmacodynamics between-subject variabilities of Ep administration in critically ill children. This 
approach ultimately leads to a valuable Ep dosing simulation which should help clinicians to determine an appropriate 
a priori dosing regimen. 



Introduction 

Inotropic agents are commonly administered to prevent 
postoperative low cardiac output syndrome (LCOS) fol- 
lowing cardiopulmonary bypass (CPB) in children under- 
going open heart surgical repair [1], According to the 
PRIMACORP study, milrinone is the first-choice drug [2]. 
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However as described in the European survey EuLoCOS- 
Paed, preventive drug therapy is highly variable. For in- 
stance, epinephrine (Ep), which is cheaper than other 
commonly used catecholamines, is also used, although 
evidenced-based data are currently lacking [3,4]. 

The amplitude of the hemodynamic response to Ep is 
difficult to predict given the multitude of factors involved 
and clinical experience suggests broad between-subject 
variability. This hemodynamic response is primarily 
dependent on Ep concentrations. However Ep pharma- 
cokinetics has been poorly evaluated in children. Fisher 
et al. suggested linear pharmacokinetics with a lower 
clearance than that reported in healthy adults, although 
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only six children were included in their study and neither 
inter-patient variability nor pharmacodynamic effects were 
described [5]. A recent adult study using population 
pharmacokinetic modeling highlighted the influence of 
bodyweight (BW) and disease severity on Ep clearance 
confirming this variability [6]. These between-subject 
disparities may be even more pronounced in children. 

Pediatric dosages of Ep are usually extrapolated from 
adult studies. However, postnatal development of cardiac 
contractility is associated with major changes in the 
modulatory effect of (3-adrenoreceptor signaling. More- 
over, differential maturation of the transduction pathways 
within the cardiomyocyte contributes to age-dependent 
changes in cardiac responsiveness and sensitivity to ago- 
nists [7]. Although much is known regarding the adult 
physiological and pharmacological effects of Ep, there are 
very few pediatric studies on Ep pharmacodynamics. Ef- 
fects of Ep infusion in children have only been described 
in the neonate, mostly in low birth-weight infants, where 
effects on heart rate (HR), mean arterial pressure (MAP), 
plasma glucose and lactate levels were observed [8]. 

The purpose of the present study was to investigate, 
using a population approach, the pharmacokinetics and 
pharmacodynamics of Ep including hemodynamic (HR, 
MAP) and metabolic effects (plasma glucose and lactate 
levels) in critically ill children undergoing surgical repair 
for congenital heart defect, following CPB, as well as as- 
sociated variability factors [9]. The effects of develop- 
mental and other factors on Ep pharmacokinetics and 
pharmacodynamics were investigated in order to better 
explain the observed between-subject variabilities and to 
ultimately suggest individualized dosage regimens. 

Materials and methods 

Setting 

This prospective study was conducted in a 14-bed surgical 
pediatric cardiovascular intensive care unit (pCVICU) of a 
tertiary teaching hospital Necker Enfants Malades, Paris 
in France from July 2011 to December 2011. The Ethics 
committee of the Necker Enfants Malades Hospital ap- 
proved the study provided that written and appropriate 
consent was obtained from the child's parent(s) after they 
were informed of the objectives. We confirm that we have 
all necessary and appropriate consent from each child's 
parents involved in the study, including consent to partici- 
pate in the study and consent to publish. 

All consecutive children aged less than 18 years, weigh- 
ing more than 1,200 g, and requiring Ep infusion following 
CPB for open heart surgery were included. Non-inclusion 
criteria were unknown initial time infusion of Ep, un- 
known time of Ep flow-rate changes or unknown time of 
blood sampling. Children were enrolled prior to the onset 
of infusion and for a period lasting 6 hours after the start 
of Ep administration. 



Intervention 

In the operating room, all children underwent endo- 
tracheal intubation and were mechanically ventilated 
under sedation, opioid treatment (midazolam and sufen- 
tanil) and neuromuscular blocking agent. Standard mon- 
itoring was used, comprising a radial or femoral artery 
catheter for measurement of systemic arterial blood 
pressure and intermittent blood sampling, a triple-lumen 
right internal jugular or femoral central venous catheter 
(CVC), and urinary bladder or rectal temperature probes. 
Normothermic CPB with intermittent warm blood cardio- 
plegia was performed in every patient during the study 
period, except in cases where deep hypothermic circula- 
tory arrest was indicated. Conventional ultrafiltration was 
performed during the CPB. 

Ep infusion was initiated in the operating room (defined 
as time = 0 minutes), in association with milrinone at the 
end of the CPB, according to the local protocol and the 
risk of developing an LCOS: risk adjustment for congeni- 
tal heart surgery 1 (RACHS-1) category, aortic cross- 
clamping duration, preoperative left ventricle dilatation, 
preoperative or intraoperative arterial pulmonary hyper- 
tension defined by intra cardiac right to left shunt or pul- 
monary arterial pressure over 2/3 mean systemic arterial 
pressure, hemodynamic instabilities (defined by a variation 
greater than 20% of HR and/or MAP) and physiological 
status [2,10]. Cases involving sepsis or preoperative myo- 
cardial dysfunction requiring inotropic support were 
excluded. 

Ep (adrenalin 1 mg-mL' 1 , Renaudin™ diluted to 
100 ug mL' 1 or 50 ug mL" 1 in Glucose 5% Baxter™, UK) 
was infused using a programmable electric syringe pump 
(DPS, Fresenius Vial™, Brezins, France) through a triple- 
lumen central venous catheter (Arrow, Teleflex , PA 
19605, USA) with a flow rate varying from 0.3 mLh 1 
to 1 mL h" 1 . The latter was adjusted for age and 
hemodynamic objectives, namely: normal HR, normal 
MAP, normal time capillary refill, normal pulse, nor- 
mal urine output (>2 mL kg ^h" 1 ), ScV02 >70% when 
measured, normal transthoracic echocardiography of 
left ventricular ejection fraction (60 to 80%) and nor- 
mal plasma lactate level (<2.2 mmol-L' 1 ) [2,11]. De- 
pending on the congenital heart defects and the repair, 
the preload was normalized based on CVP and/or left atrial 
pressure (LAP). After normalized preload, intravenous mil- 
rinone at a dose ranging from 0.3 to 0.5 ug kg^ min" 1 was 
systematically combined with Ep infusion without loading 
bolus at initiation. 

Upon arrival to the ICU, medications were adjusted by 
the bedside nurse under the direction of the medical 
team: blood transfusion to reach a hemoglobin level 
above 10 g dL" 1 , furosemide to maintain water balance 
and urine output over 2 mL-kg' h' . Adequate analgesia 
and sedation were ensured by, respectively, continuous 
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intravenous morphine or sufentanil and midazolam, 
mechanical ventilation with adequate pressure levels 
and oxygen inspired fraction and inhaled nitric oxide 
in case of pulmonary arterial hypertension. 

The daily amount of intravenous glucose was adjusted 
for age: birth to 12 months: 4 g-kg'^day' 1 , 12 months to 
48 months: 3 g-kg'^day' 1 , 48 months to 72 months: 
2.5 g kg' 1 day" 1 and over 72 months: 2 g -kg" 1 -day" 1 . 

LCOS was defined if Ep and/or milrinone were needed 
over 48 hours to maintain normal hemodynamic param- 
eters (normal HR, normal urine output, normal MAP, 
normal capillary refill time, warm extremities) without 
metabolic acidosis (standard plasma bicarbonate (HCO3) 
level less than 22 mmol-L' 1 or an increase in plasma lac- 
tate level greater than 2.2 mmol-L' 1 ) [2,11]. In this study, 
no other catecholamines or corticosteroid was used in the 
first 6 hours following open heart surgery. 

Blood sampling 

An initial blood sample (Co) was collected prior to CPB 
after which Ep infusion was initiated. A second blood 
sample (Ci) was drawn at least 60 minutes after initiat- 
ing Ep infusion. A last blood sample (C2) was drawn 
40 minutes after a change in rate flow or prior to 6 hours 
after beginning Ep infusion in the case of a constant flow 
rate. 

The 60-minute steady-state interval was chosen ac- 
cording to at least five times the Ep plasma half-life in 
healthy subjects (approximately 20 minutes) and the 
dead volume of the CVC used to infuse Ep at 0.3 to 
1 mL-h' 1 rate flow (approximately 40 minutes) [4]. 

Co was used to assess plasma levels of endogenous Ep. 
Only C 0 and Q were drawn in patients who weighed less 
than 2,500 g, according to the percentage of blood volume 
permitted by the Ethics Committee of our institution. 

Sample handling 

Blood assigned to catecholamine assays was sampled in 
EDTA-tubes placed in an ice bucket then immediately 
centrifuged at 3,000 g for 5 minutes. The plasma samples 
were then separated and immediately stored at -20°C and 
thereafter at -80°C before 24 hours running. 

Assay 

Ep concentrations were blindly determined by means of 
HPLC with colorimetric detection [12]. After thawing, 
the volume of each sample was adjusted to 4 mL by 
adding distilled water and the internal standard, dihy- 
droxybenzylamine. A 20-uL aliquot at 10°C was then 
injected into the chromatographic system comprised of 
a column (25 cm x 4.6 mm inner diameter, 5 urn 
Supelcosil LC-18 Supelco ), an electrochemical ESA 
colorimetric detector (Model Coulochem III, Eurosep™), 
dual analytic cells (ESA cell Model 5011) set at -0.05 V 



for the first detector and -0.3 V for the second detector, 
and a conditioning cell set at +0.3 V. The mobile phase, at 
1.2 mL-min' 1 , consisted of a mixture of an aqueous acidic 
buffer containing heptane sulfonic acid and acetonitrile 
(93:7 v/v). Ep calibration curves were prepared according 
to the same procedure (2.5 ug to 75 ug/4 mL distilled 
water). The measured Ep concentration in pmolmL' 1 was 
converted to ug-L' 1 . The detection threshold (defined 
by variability <10% between measurements) for HPLC was 
0.2 pmol mL' . Endogenous and exogenous Ep were 
strictly identical with regard to chromatographic detection. 

Patient data 

Baseline patient characteristics were recorded, including 
non-cardiac medical history, gender, age, BW, RACHS-1 
category, type of congenital heart defect, preoperative 
cyanotic status and left ventricular ejection fraction, dur- 
ation of CBP and aortic cross-clamping, duration of 
pCVICU stay, mechanical endotracheal ventilation dur- 
ation and death during pCVICU stay. Duration of both 
Ep and milrinone infusion were recorded. Variation of 
infused doses was recorded during the first 6 hours. 

HR (beats min' 1 ) and invasive MAP (mmHg) data were 
recorded prior to CPB, at initiation of Ep, and then every 
10 minutes for the first hour and thereafter every hour or 
less if needed during the subsequent 6 hours. Left ven- 
tricular ejection fraction (%) was measured at least once 
during the 6 hours. CVP (mmHg) was systematically 
recorded as well as LAP (mmHg) when measured. 
Temperature (°C) and urine outputs (mL-kg^-h' 1 ) were 
recorded during 6 hours following CBP. 

Plasma lactate and glucose levels (mmol-L' ) were re- 
corded before surgery and at least once thereafter during 
the following 6 hours. Arterial pH, ionized plasma calcium 
levels (mmol-L 1 ) and plasma HCO3 levels (mmol-L 1 ) 
were recorded during the first 6 hours. 

Results are expressed as raw numbers (%) or medians 
(ranges). The non-parametric Wilcoxon test was per- 
formed to compare pharmacokinetic and pharmacody- 
namic values before and under Ep infusion. P <0.05 was 
considered statistically significant. 

Pharmacokinetic-pharmacodynamic modeling 

Ep concentration time-courses were described by a one- 
compartment open model with first-order elimination 
with the parameters of elimination clearance (CL) and 
volume of distribution (V). The differential equation 
connected to this model is thus, 

dA(t)/dt = Rate - CL.C(t) (1) 

withC(t) = A(t)/V (2) 

where A(t) and C(t) denote the amount of drug and con- 
centration of drug in the body at time t. The endogenous 
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production rate, qO, was taken into account in the model 
as follows, 

C(t) = A(t)/V + qO/CL (3) 

The effect of BW was investigated in the pharmacoki- 
netic model via an allometric relationship [13]. 

P = P TYP (BW) PWR (4) 

where P, Ptyp and PWR are the individual parameter, 
typical parameter and power exponent, respectively. The 
PWR exponent was estimated in a first attempt and then 
eventually fixed to % for CL and qO terms according to 
the typical weight-based allometric rule. 

The circulating volume, Varc (L)> was related to BW 
as follows [14]. 

V ci rc = 0.08.BW (5) 

As kinetics data were best described by a one-compartment 
model with first-order elimination, the half-life (Ty a ) 
was related to V and CL as = ln2 .V/CL = 0.69.V/CL. 

The HR response, HR(t), was related to the Ep con- 
centration via an Emax model. 

HR(t) = HR 0 + (HR max -HR 0 ).C(t)/{C(t) + C 50 HR} 

(6) 

where HR max and HR 0 are respectively the maximal and 
basal HR values and C50HR the concentration that in- 
duces 50% of the maximal effect on HR. 
The MAP(t) is then expressed as: 

MAP(t) = HR(t).SV.SVR(t) + CVP (7) 

where SV, SVR and CVP represent stroke volume, sys- 
temic vascular resistance and central venous pressure, re- 
spectively [15]. As SV and SVR were not known, the SV. 
SVR product variation was estimated via the function. 

SV.SVR(t) = SV.SVRo + (SV.SVR max -SV.SVR 0 ).C(t)/ 
{C(t) + QoSV.SVR} 

(8) 

where SV SVR 0 , SVSVR max and C 50 SV SVR respectively 
denote the SVSVR product basal value, the product's 
maximal value and the concentration that induces 50% 
of the maximal effect on SV SVR. 

Plasma glucose and lactate, G(t) and L(t), variations 
were modeled by a turnover model in which the stimula- 
tion of plasma glucose production, S(t), was related to 
Ep concentration as follows. 

dG(t)/dt = Rgly ■ S ( t) — 1<gly • G ( t) (9) 
with 

S(t) = 1 + G max C(t)/{C(t) + C 50 GLY} (10) 



where R GL y and 1< GL y represent the plasma glucose zero- 
order rate production and first-order elimination rate 
constant. G max and C 50 GLY denote respectively the 
maximal stimulation effect and the Ep concentration 
that produces 50% of the maximal response. 

The rate of change in plasma lactate level, dL(t)/dt, 
was related to the plasma glucose level variation rate as: 

dL(t)/dt = k GLY .G(t)-k LA c.L(t) (11) 

where 1<lac is the plasma lactate elimination rate constant. 

Before Ep infusion, the systems are assumed to be at 
steady-state, 

G(0) = GLY 0 , L(0) = LAC 0 , then 

1<qly an d 1<lac are 

k GLY = Rgly/GLYq (12) 

1<LAC = kGLY- GLY0/LAC0 (13) 

where GLY 0 and LAC 0 denote respectively, basal plasma 
glucose and lactate levels. 

Population pharmacokinetic-pharmacodynamic analysis 

Drug concentrations and responses were analyzed using a 
population approach, that is, a non-linear mixed-effect mod- 
eling approach. Data were analyzed using the MONOLIX 
software version 4.13 s [16] and the SAEM algorithm [17]. 
Differential equations were written in an MLXTRAN script 
file in MONOLIX to estimate the parameters. Residual 
variabilities were described by additive, proportional or ex- 
ponential error models depending on the observation. An 
exponential model was used for between-subject variability 
(BSV). The effect of a covariate on a structural parameter 
was retained if it caused a decrease in the Bayesian informa- 
tion criterion (BIC) and/or reduced the corresponding BSV 
with P <0.05. Only covariates with a plausible effect on phar- 
macokinetic and pharmacodynamic parameters were inves- 
tigated. The main covariates of interest in this pediatric 
population were BW and age. 

Visual predictive check (VPC) evaluation 

Plasma Ep concentration, HR, MAP, plasma glucose and 
lactate-level time course was simulated from the respect- 
ive final population model and compared with the ob- 
served data to evaluate the predictive performance of 
the model. The vector of pharmacokinetic parameters 
from 400 replicates of the database was simulated using 
the final model. Each vector parameter was drawn in a 
log-normal distribution with a variance corresponding to 
the previously estimated BSV. A simulated residual error 
was added to each simulated concentration. The 5th, 50th 
and 95th percentiles of the simulated dependent variables 
at each time point were then overlaid on the observed 
data and a visual inspection was performed. Because the 
patients received different Ep regimens, the Uppsala 
correction was used to produce the VPC plots [18]. 
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Evaluation and validation 

Diagnostic graphics were used for evaluation of the 
goodness-of-fit. Concentration and effects profiles were 
simulated and compared with the observed data with 
the aid of the VPC in order to validate the model. 

Results 

Patients 

A total of 55 children were initially enrolled, of which 16 
patients were subsequently excluded: 6 because of incom- 
plete parental consent, 7 because of missing Q and C2 
blood samples and 2 because of hemolysis. 

Hence, 39 children were included in the study. Co sam- 
ples were obtained in 33 patients, Q in all children and C2 
in 25 children for a total of 97 observations. Hemodynamic 
data (HR, MAP) and metabolic effects of Ep infusion 
(plasma lactate and glucose levels) were available in 38 chil- 
dren with 434, 464, 101 and 140 observations, respectively. 

Five premature children with a gestational age <37 weeks 
(n = 1 at 33 weeks, n = 1 at 34 weeks, and n = 3 at 36 weeks) 
were recorded. Chromosomal disorders were reported in 
eight children (n = 3 with Down syndrome, n = 1 with 
di George syndrome, n = 1 with Loeys-Dietz syndrome, 
and n = 1 with Noonan syndrome n = 2 with suspected 
Noonan syndrome). Respiratory disorders were noted 
in seven patients (n = 3 with asthma, n = 3 with laryn- 
gotracheomalacia, and n = 1 with chronic aspiration 
pneumonia) and malnutrition (<2 SD) was observed in 
nineteen children. 

Six children were treated before open heart surgery 
with converting enzyme inhibitors because of left ven- 
tricular dilatation, seven were treated with prostaglan- 
dins because of ductus arteriosus-dependent heart disease, 
and |3-blockers were co-administered to three children be- 
cause of obstruction of the left ventricular outflow track. 

All children had transthoracic echocardiography prior 
to CPB; left ventricular ejection fraction was evaluated at 
60% (25 to 78) with normal values (>50%) observed in 
34 patients (87%). Ventricular diastolic function was not 
assessed. Eleven children were cyanotic (Sp02 <90%) prior 
to the surgery because of their congenital heart disease. 
Open heart surgeries were as follows: arterial switch (n = 9), 
repair of complete atrioventricular canal (n = 6), repair of 
ventricular septal defect (n = 6), total repair of tetralogy of 
Fallot (n = 4), pulmonary valvuloplasty (n = 3), repair of co- 
arctation and ventricular septal defect closure (n = 2), repair 
of interrupted aortic arch with ventricular septal defect 
closure (n = 2), repair of double-outlet right ventricle 
(n = 2), repair of pulmonary artery stenosis (n = 2), pul- 
monary valvular replacement (n = 2), and repair of truncus 
arteriosus (n = 1). Surgery with deep hypothermic circula- 
tory arrest was necessary in 1 1 children. 

In the operating room all children required the follow- 
ing: a red blood cell transfusion, fresh-frozen plasma 



administration, neuromuscular blocking agents, and 
hypnotic and opioid drug infusion. Ultrafiltration of 
650 mL (250 to 1,200) of fluid during CPB was per- 
formed to achieve a negative fluid balance and hematocrit 
at 44% (35 to 47). Milrinone and Ep were initiated just be- 
fore weaning from CPB, with an infused Ep dose of 
0.07 ugkg^min 1 (0.01 to 0.23) and an infusion milrinone 
dose of 0.5 |ig-kg" •min" (0.2 to 0.7). Milrinone infusion 
was stopped after 1.5 days (1 to 13). 

Delayed sternal closure occurred in four patients. Post- 
operative left ventricular ejection fraction under Ep and 
milrinone infusion was at 60% (30 to 70) with normal 
values (>50%) observed in 34 patients (87%). CVP and 
LAP were 11 mmHg (8 to 15) and 8 mmHg (6 to 14), 
respectively. Four children exhibited supraventricular 
tachycardia, one had ventricular tachycardia and five had 
a transient atrioventricular block, which required exter- 
nal cardiac pacing. Urine output was 4.5 mL-kg h' 1 (0.8 
to 7.5) and pH, plasma HCO3 (mmol-L 1 ) and plasma 
ionized calcium (mmol L 1 ) levels were: 7.39 (7.27 to 
7.45); 24 (20 to 26) and 1.28 (1.12 to 1.5), respectively. 
All children required diuretics, whereas none were under 
corticosteroid therapy during the 6 hours following sur- 
gery. Ten children needed inhaled nitric oxide for pulmon- 
ary arterial hypertension during the postoperative course. 
Endotracheal mechanical ventilation was performed for all 
patients during 2.1 days (1 to 17). None of the patients re- 
quired renal replacement therapy and none had liver injury 
according to prothrombin activity and/or Factor V levels 
(lower than 50% for at least 24 hours). None of the chil- 
dren died during their pCVICU stay. Finally, nine children 
(23%) had LCOS according to the classical definition. 
Table 1 summarizes overall patient characteristics. 

Epinephrine pharmacokinetics 

The increase in Ep concentration during infusion was sig- 
nificant: 2.94 ug-L~ (0.37 to 71) compared to baseline Ep 
concentration, 0.062 ug-L" 1 (0.037 to 0.25) (P <0.001). A 
one-compartment open model with linear elimination ad- 
equately described the Ep time courses. An additional 
movie file shows this in more detail (see Additional file 1). 
The pharmacokinetic parameters were V, CL and qO. The 
residual variability was ascribed to a proportional 
model. BW was the main covariate influencing CL and 
qO (.P <0.001). Both PWR estimates for qO and CL were 
close to 1 (0.98 and 0.985), however CL and qO were 
poorly estimated (relative standard errors near 50%). 
Moreover, BIC decreased from 214.2 to 205.0 when 
PWR was estimated and further decreased to 199.0 
when PWR was fixed to %. Also, there was no visible 
or significant difference between the two models on 
the observed versus predicted plots. An additional 
movie file shows this in more detail (see Additional 
files 2 and 3). 



Oualha et al. Critical Care 2014, 18:R23 
http://ccforum.eom/content/1 8/1 /R23 



Page 6 of 14 



Table 1 Patient characteristics and baseline kinetic and 



Table 2 Population pharmacokinetic parameters 



dynamic data 




Parameters 


Estimate 


RSE (%) 


Patient characteristics (n = 39) 


Values 


6 CL (L-h^-kg- 1 ) 


2 


17 


Demographics 




e BW (CL(BWi) = e CL x bw? /4 ) 


0.75 (fixed) 


NA 


Age, months 


3.9 (0.1 to 189) 


9 q0 (Mg-h-'V) 


0.15 


19 


Gender, male, n (%) 


26 (66.6%) 


eBwtqOtBWiHGqoX BWf /4 ) 


0.75 (fixed) 


NA 


Body weight, kg 


4.5 (2.5 to 58) 


V (L) for a 10 Kg individual 


0.8 


NA 


Preoperative physiological profile 




Ty 2 (min), for a 10 Kg individual 


3 


NA 


RACHS-1, categories, median (range) 


3 (2 to 4) 


Hcl (square root of to CL ) 


1 


14 


Category 2, n 


16 


n q o (square root of to 2 q0 ) 


1.1 


13 


Category 3, n 


17 


Residual variability (proportional) 


0.3 


15 


Category 4, n 


6 


Correlation (n CL , n q o) 


0.88 


5 



Prothrombin activity (%) 80 (50 to 1 00) 

Creatinine clearance, mL- min" 1 -1.73 m (2) _1 91 (22 to 200) 
Per-operative course 

Duration of CPB, minutes 1 07 (52 to 222) 

Aortic cross-clamping time, minutes 64 (9 to 140) 
Postoperative course 

Ep infusion duration, days 1.5 (1 to 13) 

pCVICU length of stay, days 3 (2 to 23) 
Baseline kinetic and dynamic parameters 

Ep concentration, ug-L 1 0.062 (0.037 to 0.25) 

Heart rate, beats- min" 1 135 (70 to 180) 

Mean arterial pressure, mmHg 51 (25 to 65) 

Plasma glucose level, mmol-L 6.2 (4 to 10.1) 

Plasma lactate level, mmol-L 1 1 (0.5 to 3) 

RACHS, risk adjustment for congenital heart surgery; n, number; CPB, 
cardiopulmonary bypass; min, minutes; Ep, epinephrine; pCVICU, pediatric 
cardiovascular intensive care unit. Values represent median and the 
presentation of the value in parentheses defined the range from (min to max). 

Hence, these power values were fixed to % according 
to the BW-based allometric rule. V was assumed to be 
equal to the circulating volume because V could not be 
accurately estimated and because of the hydrophilic na- 
ture of Ep. No other covariate (gender, pH, temperature, 
RACHS-1 categories, creatinine clearance, gestational 
age, cyanosis, malnutrition) influenced the pharmacokin- 
etics. The final relationship for Ep CL and qO was: 



CL(BWi) = Sex x (BWif 4 and qO(BWi) = 6 q0 x (BWi )**, 
then 

e CL (L.h^.kg 1 ) = 2, e q0 (ug.h-^kg- 1 ) = o.i5 

where 8 CL and 6 q0 are typical unit clearance and en- 
dogenous production rate. For an individual weighing 
10 kg, Ep CL, qO, V and T% were: 

CL (10 kg) = 2 x 10 y * = 11.2 L-h"\ qO (10 kg) = 0.15 x 
10* = 0.84 ug h" 1 , V(10 kg) = 0.08 x 10 = 0.8 L and 
T H (10 kg) = Qn2 x V(10 kg) /CL(10 kg)) x 60 = 3 minutes. 



The volume of distribution of epinephrine was ascribed to the circulating 
volume, estimated as a function of bodyweight (see Materials and methods). 
CL, elimination clearance; qO, endogenous production rate; V, volume of 
distribution; CV, circulating volume; 6 CL , typical unit clearance; 9 q0 ,typical unit 
endogenous production rate; 6 BW , bodyweight influential parameter; Ty 2r half-life. 
RSE, relative standard error; r), between-subject variability (BSV); BW, bodyweight; 
NA, not applicable. 

Table 2 summarizes the final population estimates. All 
parameters were estimated with good precision. Figure 1 
depicts the VPC and shows that the observed concentra- 
tions were well centered around the simulated median 
predictions. 

Epinephrine hemodynamics 

After initiation of Ep infusion, (i) HR increased significantly 
from 135 beats min" 1 (70 to 180) to 159 beats min" 1 (80 to 
212) (variance (w) = 2563; P = 2.10~ 8 ), (ii) MAP values 
increased significantly from 51 mmHg (25 to 65) to 
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Figure 1 Prediction corrected-visual predictive check (PC-VPC) 
for epinephrine concentrations versus time in minutes. The 

green lines depict the 5th, 50th and 95th percentiles of observed 
data; the areas represent the 95% confidence interval around the 
simulated percentiles. Blue color represent the 5th and the 95th 
percentile of the predicted concentration versus time and pink color 
represent the median predicted concentration versus time. Time 
0 min represents the starting time of epinephrine infusion. 
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66 mmHg (30 to 94) (variance (w) = 2613; P = 5.10" 11 ) 
(Figure 2). The variations in HR and MAP as a func- 
tion of Ep concentration were well explained by the 
Emax models, expressed by equations (6) to (8). The 
residual variability was ascribed to a proportional model. 
BSVs were estimated for HR 0 , C 50 HR, SV.SVR 0 and SV. 
SVR max . Age was the main covariate influencing HR 0 and 
SV-SVRo (P <0.001) where HR oi = HR 0 x age; 00612 and 
SV.SVRoi = SV-SVRo x age 0094 respectively. Including age 
in the model dramatically decreased the BIC (from 5,998 



to 5,971) and improved the curve-fitting of the model. In 
addition, RACHS-1 category was significant (P = 0.04) in 
the estimation of SV-SVR max : 0.44 and 0.26 for RACHS-1 
categories 2 and 3 to 4 respectively; the BSVs for HR 0 , 
C 50 HR, SV-SVRo and SV SVR max varied from 0.19, 1.0, 
0.25 and 0.32 (basic model) to 0.14, 1.22, 0.13 and 0.23 
(final model). Also, the BIC decreased from 5,971 (includ- 
ing age) to 5,965 (final model). 

No other covariate improved the model (including dur- 
ation of aortic cross-clamping, duration of CPB, deep 



B 



during infusion at steady-state 



during infusion at steady-state 



i i 1 1 

basal level during infusion at steady-state basal level during infusion at steady-state 

Figure 2 Box and whisker plots of hemodynamic and metabolic data before and during epinephrine infusion. Box and whisker plots of 
(i) heart rate, beats- min" 1 (A) before epinephrine infusion (minimum = 70; maximum (max) = 180; median = 135) and during epinephrine infusion 
(minimum = 70; max = 212; median = 159) (variance (w) = 2,563; P= 2.10" 8 ); (ii) mean arterial pressure, mmHg (B) before epinephrine infusion 
(minimum = 25; max = 65; median = 51) and during epinephrine infusion (minimum = 30; max = 94; median = 66) (variance (w) =2613; P = 5.10" 11 ); 
(iii) plasma glucose, mmol-L" 1 (C) before epinephrine infusion (minimum = 4; max= 10.1; median = 6.2) and during epinephrine infusion 
(minimum = 4.9; max = 23.3; median = 9.85) (variance (w) = 339; P = 6.10" 9 ) and (iv) lactate levels, mmol-L" 1 (D) before epinephrine infusion 
(minimum = 0.5; max = 3; median = 1) and during epinephrine infusion (minimum = 0.9; max = 7.1; median = 2.1) (variance (w) = 218; P = 3 . 1 0" 1 °) . 
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hypothermia and milrinone dosage). The final population 
parameters are summarized in Table 3. The VPC plots in 
Figure 3 show that the observed HR and MAP values are 
well-centered around the predicted median of the model. 

Metabolic effects of epinephrine 

Both plasma glucose and lactate levels increased signifi- 
cantly after the initiation of Ep infusion from 6.2 mmol-I/ 1 
(4.0 to 10.1) and 1 mmol.L" 1 (0.5 to 3) to 9.85 mrn.ol.L~ 1 (4.9 



Table 3 Hemodynamic and metabolic population 
parameters 



Hemodynamic population parameters 


Estimate 


RSE (%) 


HR 0 (b-min"') 


133 


3 


e a9 e (HR oi = HR 0 xageL a06 ') 


-0.061 


18 


HR max (b-min" 1 ) 


180 


3 


C 50 HR (ug-L" 1 ) 


5.71 


37 


SV-SVR 0 


0.31 


3 


9 age (SV SVRoi = SV SVR 0 x age?' 094 ) 


0.094 


14 


SV-SVR max , RACHS-1 = category 2 


0.44 


20 


SV.SVR max , RACHS-1 = categories 3,4 


0.26 


20 


C 50 SVSVR (ug.L 1 ) 


18 


59 


Hhro (square root of to 2 HR 0 ) 


0.14 


14 


Hcso hr (square root of to 2 c 50 HR ) 


1.22 


17 


Hsv.svro (square root of id 2 sv.svro) 


0.13 


17 


HsV.SVRmax (square root of (J 2 SV.SVRmax) 


0.23 


42 


Residual variability (proportional) 






HR 


0.08 


4 


MAP 


0.16 


4 


Metabolic population parameters 






GLY 0 (mmol-L 1 ) 


5.46 


5 


("max 


1.69 


6 


R G ly (mmol-L" 1 .min -1 ) 


0.04 


25 


C 50 GLY (ug-L 1 ) 


0.52 


9 


LAC 0 (mmol-L 1 ) 


1.23 


/ 


Hglyo (square root of co 2 GLY 0 ) 


0.21 


23 


HGmax (square root of io 2 Gmax ) 


0.213 


26 


n RG iy (square root of cu 2 RGLY ) 


1 


21 


Hlaco (square root of uj 2 LAC0 ) 


0.33 


18 


Residual variability (constant additive) 






GLY (mmol-L 1 ) 


2.23 


5 


LAC (mmol-L 1 ) 


0.5 


ll 



HR 0 , basal heart rate; HR max , maximal HR; C 50 HR, Ep concentration producing 
50% of HR max ; SV SVR, product of stoke volume by systemic vascular 
resistances; SVSVR 0 , basal SV SVR value; SV.SVR max , maximal SV SVR value; 
C 50 SV-SVR Ep concentration producing 50% of SVSVR max ; RACHS-1, risk ad- 
justment for congenital heart surgery 1; GLY 0 , basal plasma glucose level; 
R GLY , glucose zero-order production rate; G max , relative maximal increase in 
R GLY ; C 50 GLY, Ep concentration that produces 50% of the maximal response 
on plasma glucose level; LAQ, basal plasma lactate level; RSE, relative standard 
error; q, between-subject variability (BSV); 6 age , age influential parameter. 



to 23.3) [variance (w) = 339; P=6.10" 9 ] and 2.1 mmol-L" 1 
(0.9 to 7.1) (variance (w) = 218; P = 3.10" 10 ), respectively 
(Figure 2). Assuming that Ep stimulated the glucose pro- 
duction rate, the turnover models expressed by equations 
(9) to (13) effectively described the variations in plasma 
glucose (G(t)) and lactate levels (L(t)). There was no signifi- 
cant covariate effect, including those of exogenous glucose 
supply, age or BW. The residual variability was ascribed to 
a constant additive model. BSVs were estimated for GLY 0 , 
Rgly, G max and LACn. All parameters were well-estimated 
with low relative standard errors. Table 3 summarizes the 
population estimates. The VPC plots in Figure 4 show that 
the observed plasma glucose and lactate levels are well- 
centered around the predicted median of the model. 

Epinephrine dosing simulations 

Using the hemodynamic model, the effects of various infu- 
sion rates of Ep on HR and MAP were assessed as a func- 
tion of age and BW. As shown in Figure 5, the increase in 
Ep concentration versus infusion rate was linear although 
the increases in HR and MAP were curvilinear, due to 
the Ep concentration-Emax model. Similarly, Figure 6 
shows the metabolic responses for a child weighing 5 kg 
with three infusion rates: 0.02, 0.1 and 0.25 ug-kg-min" 1 . 

Discussion 

Little is known about the pharmacokinetics and pharma- 
codynamics of Ep in children although this drug is used 
in the pediatric population. The present study using a 
population approach adequately described the kinetics, 
hemodynamic and metabolic effects of Ep in critically ill 
children, highlighting the between-subject variabilities 
which were well explained by age and BW. 

Epinephrine pharmacokinetics 

A one-compartment open model with linear elimination 
adequately described the data as previously reported 
[3,5,19]. The effect of BW using the allometric scale on 
clearance and Ep endogenous production improved the 
model and partly explained the between-subject variabil- 
ity [13]. This was not unexpected since endogenous rates 
of production and clearance of Ep are dependent on en- 
zymatic maturation, both of which are related to age and 
BW [20]. With regard to endogenous Ep rate produc- 
tion, given that the concentrations observed following 
the infusion were well above 10-fold the baseline con- 
centrations (approximately 50-fold), the contribution of 
possible variations in endogenous production was as- 
sumed to be negligible during the infusion. We could 
not adequately estimate volume of distribution because 
Ep concentration was measured only in the steady state; 
however, adjusting the volume of distribution to the 
circulating volume is justified considering the hydro- 
philic nature of Ep. We did not find any effect of 
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Figure 3 Prediction corrected-visual predictive check (PC-VPC) for heart rate (A) and mean arterial pressure (B) observations versus 
time in minutes. The green lines depict the 5th, 50th and 95th percentiles of observed data; the areas represent the 95% confidence interval 
around the simulated percentiles. Blue color represent the 5th and the 95th percentile of the predicted HR and MAP versus time and pink color 
represent the median predicted HR and MAP versus time. Time 0 min represents the starting time of epinephrine infusion. 



creatinine clearance, as only 10% of Ep is excreted 
unchanged via the renal route and is mainly and rapidly 
inactivated by either methylation via the effect of 



catechol-O-methyl transferase or oxidative deamination 
by monoamine oxidase into inactive metabolites excreted 
by the kidney [21]. In contrast to the study of Abboud 
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Figure 4 Prediction corrected-visual predictive check (PC-VPC) for plasma glucose level (A) and plasma lactate level (B) observations 
versus time in minutes. The green lines depict the 5th, 50th and 95th percentiles of observed data; the areas represent the 95% confidence interval 
around the simulated percentiles. Blue color represent the 5th and the 95th percentile of the predicted plasma glucose level and lactate level versus time and 
pink color represent the median predicted plasma glucose level and lactate level versus time. Time 0 min represents the starting time of epinephrine infusion. 



et al., neither RACHS-1 categories nor duration of CBP 
or aortic cross-clamping, which reflect the severity of 
illness, were found to be significant, possibly because of 
the small sample size and the difference in patient age 
and illness groups between the studies [6]. 



Epinephrine - hemodynamic effects 

To the best of our knowledge, this is the first study in 
which the hemodynamic responses to Ep in preventing 
LCOS were modeled. Only HR and MAP were recorded 
in this study. In adult volunteers, as in critically ill 
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Figure 5 Dosing simulations depicting the increases in epinephrine concentration and hemodynamic responses as a function of 
administered infusion rate in children of different bodyweights and ages for patients with a risk adjustment for congenital heart 
surgery (RACHS)-1 category of 2. 



patients, Ep increases HR as well as MAP [8,22,23]. The 
predominant effect of Ep when administered at low dose 
is mediated by (3-adrenergic receptors, which increase 
HR and SV [23-25]. The resulting hemodynamic re- 
sponse may differ in children because of (i) the relative 
immaturity of the myocardium, which precludes a sig- 
nificant increase in SV [26] and (ii) a variation in pi and 
|32-adrenergic receptor density related to age [27]. Fur- 
thermore, newborn and infant cardiac outputs are more 
dependent on an increase in HR [25]. This is confirmed 
by our hemodynamic model, which highlighted the pri- 
mary role of HR, which in turn influences MAP, that is, 
the C 50 estimate for the SVSVR product was three 
times as high as that estimated for HR, reflecting a 
much greater sensibility of the HR response to Ep. 

The Emax models were effective in relating both HR 
and MAP responses to Ep concentration [28]. Measuring 
cardiac output is hazardous in children undergoing re- 
pair of congenital heart disease because of residual ven- 
tricular and/or auricular shunt. Hence, we could only 
estimate the SV SVR product that relates MAP to HR. 



Moreover, as expected, age was found to be a significant 
covariate that dramatically improved the model: HR was 
negatively related to age whereas the SV SVR product 
increased with age. This latter finding is easy to explain 
because both SV and MAP increase with age [29]. Like- 
wise, RACHS-1 categories 3 and 4 decreased the max- 
imal MAP response (SV SVR m£tx product). This effect 
may be related to the role of systemic inflammatory syn- 
drome following CB, which alters myocardial and vascu- 
lar response to Ep [11]. The impact of temperature and 
pH on HR was also investigated but was not found 
significant. 

As milrinone was infused in all children, this could 
have had a confounding effect on hemodynamic re- 
sponses. There are conflicting data in the literature on 
the effect of milrinone on HR, namely an increased or 
unchanged HR [30]. In the present study, none of the 
parameters of the hemodynamic model were found to be 
influenced by a dose-dependent effect of milrinone. 
Moreover, any possible confounding effect would be neg- 
ligible since an immediate increase in HR was observed 
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after Ep initiation which is not compatible with the de- 
layed response to milrinone [31]. 

Epinephrine - metabolic effects 

The turnover model herein provided a valid relationship 
between Ep concentration and plasma glucose and lac- 
tate levels, as these metabolic responses are dependent 
on stimulation of glycogenolysis via the activation of 
pVadrenergic receptors [32]. During Ep infusion, the in- 
creases in plasma glucose and lactate levels were signifi- 
cant, albeit delayed as compared to the hemodynamic 
responses. Indeed, exogenous Ep has previously been 
shown to increase plasma glucose and lactate levels 
[8,33]. Lactate is mainly produced via the anaerobic 
glycolytic breakdown of glucose to pyruvate [32]. An ex- 
cessive vasoconstriction mediated by cc-adrenergic recep- 
tors results in lactate accumulation due to limited 
oxygen supply [34]. However, lactate may also accumu- 
late during accelerated aerobic glycolysis driven by Ep 
[35], and it is unlikely that the rise in plasma lactate 



level is due to an excessive vasoconstrictor effect via 
Ep a-adrenergic receptor stimulation, given the low 
doses administered. Lastly, we believe that these in- 
creases are strongly related to Ep without confounding 
factors because (i) there was no significant effect of the 
exogenous supply of glucose, (ii) there were no other 
potential hyperglycemic treatments, such as corticoste- 
roids and (iii) milrinone does not elevate glucose and/ 
or lactate levels [36]. 

Epinephrine dosing simulations 

Using the final model, it was possible to highlight the 
differences in responses to a same infusion rate accord- 
ing to age or BW. Therefore, these simulations allow the 
determination of an a priori dosing schedule, for specific 
BW and age, to produce a suitable increase in HR and 
MAP. Interestingly, these plots clearly show that the 
amplitude of HR or MAP increase following various Ep 
infusion rates is related to the child's BW, that is, the 
lower the BW, the smaller the amplitude of increase. 
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Limitations of the study 

The small sample size likely limited the identification of 
other significant covariates that could affect either the 
pharmacokinetics or the responses to Ep. We were not 
able to evidence a hemodynamic effect of milrinone in 
the model; however, we cannot ignore the potential ef- 
fect of the latter and must assume that the modeling of 
the hemodynamic effects of Ep, including its simula- 
tions, implicitly take into account the effects of milri- 
none. Furthermore, such simulations need to be confirmed 
in a future clinical study. Finally, as only children who 
underwent an open heart surgery with CPB were included, 
our results cannot be easily extended to patients with other 
circulatory failure etiologies. 

Conclusion 

This original study on pharmacokinetics, hemodynamic 
and metabolic effects of Ep to prevent postoperative 
LCOS in children highlights as expected, clear between- 
subject variability related to the substantial role of age 
and BW. Taking into account these individual character- 
istics should help clinicians in determining an appropri- 
ate a priori dosing regimen. 

Key messages 

• In critically ill cardiac postoperative children, lower 
bodyweight was associated with lower epinephrine 
clearance. 

• Differences of hemodynamic responses to 
epinephrine were related to age and bodyweight: the 
lower the bodyweight, the smaller the amplitude of 
heart rate and mean arterial pressure increase. 

• Increase of plasma glucose and lactate levels was 
related to epinephrine concentration without any 
effect of age, bodyweight or exogenous glucose 
supply. 

• Epinephrine dosing simulations should help the 
clinician in determining an appropriate a priori 
dosing regimen. 

Additional files 
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A: amount; BIC: Bayesian information criterion; BSV: between-subject 
variability; BW: bodyweight; C: concentration; C 50 GLY: epinephrine 



concentration that produces 50% of the maximal response on plasma 
glucose level; C 50 HR: epinephrine concentration that induces 50% of the 
maximal effect on heart rate; C 5 oSV'SVR: epinephrine concentration that 
induces 50% of the maximal effect on stroke volume-systemic vascular 
resistance; Ore: circulating volume; CL: clearance; CPB: cardiopulmonary 
bypass; CVC: central venous catheter; CVP: central venous pressure; 
Ep: epinephrine; GLY: plasma glucose level; G max : maximal increase on 
glucose zero-order rate production; HPLC: high-performance liquid 
chromatography; HR: heart rate; HR 0 : basal heart rate; HR max : maximal 
heart rate value; k GLY : first-order elimination rates; k LAC : plasma lactate 
elimination rate; LAC: plasma lactate level; LCOS: low cardiac output 
syndrome; MAP: mean arterial pressure; pCVICU: pediatric cardiovascular 
intensive care unit; PC-VPC: prediction corrected-visual predict check; 
PWR: power exponent; qO: endogenous production rate of epinephrine; 
RACHS-1 : risk adjustment for congenital heart surgery 1; R G ly: plasma glucose 
zero-order rate production; SD: standard deviation; SV: stroke volume; 
SV SVR 0 : basal value of stroke volume-systemic vascular resistance product; 
SVSVR max : maximal value of stroke volume-systemic vascular resistance product; 
SVR: systemic vascular resistance; t: time; T /2 : half-life; V: volume of distribution; 
6 CL : typical unit clearance; 8 q0 : typical unit endogenous production rate. 
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